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Inelastic  vector  soliton  collisions:  a  lattice-based 
quantum  representation 
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Lattice-based  quantum  algorithms  are  developed  for  vector  soliton  collisions  in  the 
completely  integrable  Manakov  equations,  a  system  of  coupled  nonlinear  Schrodinger 
(coupled-NLS)  equations  that  describe  the  propagation  of  pulses  in  a  birefringent 
fibre  of  unity  cross-phase  modulation  factor.  Under  appropriate  conditions  the  exact 
2-soliton  vector  solutions  yield  inelastic  soliton  collisions,  in  agreement  with  the 
theoretical  predictions  of  Radhakrishnan  et  al.  (1997  Phys.  Rev.  E56,  2213).  For 
linearly  birefringent  fibres,  quasi-elastic  solitary-wave  collisions  are  obtained  with 
emission  of  radiation.  In  a  coupled  integrable  turbulent  NLS  system,  soliton  turbu¬ 
lence  is  found  with  mode  intensity  spectrum  scaling  as  k~6. 

Keywords:  vector  solitons;  Manakov  equations;  inelastic  soliton  collisions;  quantum 
.  lattice  gas  algorithms;  turbulence 


1.  Introduction 

There  is  much  to  be  gained  in  using  micro  or  meso  descriptions  for  macroscopic 
processes.  For  example,  the  straightforward  solution  of  the  nonlinear  conservation 
equations  (e.g.  the  Navier-Stokes  equation  for  fluid  turbulence)  is  faced  with  the 
problem  of  determining  an  accurate  resolution  of  nonlinear  convective  derivatives. 
This  so-called  Riemann  problem  is  computationally  very  expensive  (Peyret  1996). 
However,  if  one  proceeds  to  a  mesoscopic  description  using  the  lattice  Boltzmann 
equation  with  the  linearized  BGK  collision  operator  (Succi  2001),  then  one  finishes 
with  a  simple  highly  parallelized  algorithm  which  in  the  long-wavelength  long-time 
limit  recovers  the  original  nonlinear  macroscopic  equations.  The  nonlinear  convective 
derivatives  of  the  macroscopic  system  are  replaced  by  simple  linear  (kinetic)  advec- 
tion  in  the  mesoscopic  lattice  Boltzmann  scheme.  Alternatively,  one  could  consider 
a  microscopic  representation  such  as  a  lattice  gas  model  with  its  detailed  kinetic 
two-body  and  three-body  lattice  collision  rules  (Frisch  et  al.  1987). 

The  nonlinear  Schrodinger  (NLS)  equation  has  played  a  pivotal  role  in  soliton 
physics  and,  in  particular,  optical  solitons  (Kivshar  &  Agrawal  2003).  Here  we 
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develop  at  the  microscopic  scales  a  quantum  lattice  gas  (QLG)  algorithm  to  exam¬ 
ine  vector  solitons  in  certain  classes  of  NLS  equations.  The  impetus  behind  quan¬ 
tum  algorithms  lies  in  their  exponential  speed-up  over  classical  algorithms  (as.  for 
example,  illustrated  by  the  Shor  (1994)  quantum  factoring  algorithm  over  classical 
factoring  schemes).  The  fundamental  change  in  moving  from  a  classical  to  quantum 
algorithm  is  to  replace  the  classical  binary  bits  ‘0’  and  T  by  the  quantum  states 
•|0)  and  |1).  It  is  the  entanglement  of  multiple  qubits,  a  purely  quantum  property 
that  has  no  counterpart  in  classical  physics,  that  can  lead  to  exponential  speed¬ 
up  (Nielsen  &  Chuang  2000).  In  the  quantum  algorithms  presented  here,  this  qubit 
entanglement  is  achieved  by  a  unitary  collision  operator  coupling  on-site  qubits.  This 
entanglement  is  then  spread  throughout  the  lattice  by  unitary  streaming  operators. 
Quantum  coherence  needs  to  be  retained  during  the  application  of  a  sequence  of 
these  collide-stream  unitary  operators,  until  quantum  measurement  is  required  to 
determine  the  amplitude  of  the  wave  function. 

Soliton  transmission  in  optical  fibres  is  an  important  aspect  of  long-distance 
(transoceanic)  communication.  However,  optical  fibres  are  birefringent,  and  a  single¬ 
mode  fibre  can  support  two  distinct  orthogonal  polarizations:  the  O-mode,  which 
has  a  constant  refractive  index  along  the  ray  path,  and  the  X-mode,  whose  refractive 
index  varies  along  the  ray  path.  Hence  pulses  with  components  in  these  orthogonal 
polarizations  will  travel  at  slightly  different  speeds.  It  can  be  shown  (Lakshmanan  & 
Kanna  2001)  that  the  slowly  varying  amplitudes  of  these  modes  satisfy  the  coupled- 
NLS  equations 

i +  2/4|<3i|2  +B\Q2\2}  ■  Qi  =  0, 

2  (1-1) 

+  +  +  ’  #2  =  0, 

where  p  is  a  positive  coefficient, 


2  +  2  sin2  9 
2  +  cos2  9 


is  the  cross-phase  modulation  coefficient  and  9  is  the  birefringence  ellipticity  angle 
(Menyuk  1989).  In  general,  the  coupled-NLS  system  equation  (1.1)  is  non-integrable. 
However,  for  the  special  case  B  =  1,  equation  (1.1)  is  completely  integrable  and  is 
known  as  the  Manakov  (1974)  system.  The  Manakov  equations  also  describe  the 
effects  of  mean  random  birefringence  on  the  propagation  of  a  circularly  polarized 
pulse  down  a  real  fibre. 

In  §  2,  we  shall  briefly  review  the  exact  analytic  solutions  for  1-soliton  and  2-solitons 
in  the  scalar  and  vector  NLS  equations.  While  the  scalar  NLS  soliton-collision  proper¬ 
ties  are  well  known — the  solitons  retain  their  exact  shape  and  speed  after  the  collision 
and  suffer  only  a  small  spatial  shift  as  the  signature  that  a  collision  has  occurred — it 
has  only  recently  been  discovered  (Radhakrishnan  et  al.  1997)  that  vector  2-soliton 
collisions  for  the  fully  integrable  Manakov  system  (equation  (1.1)  with  B  —  1)  can 
suffer  inelastic  collisions.  This  has  led  to  a  flurry  of  activity  (Jakubowski  et  al.  1998; 
Steiglitz  2001a,  b)  in  designing  sequences  of  such  soliton  collisions  that  can  affect  logic 
operations  (e.g.  controlled-NOT  gates)  and  so  manipulate  solitons  to  perform  com¬ 
putations.  In  §  3,  we  present  the  QLG  representation  for  the  coupled-NLS  system. 
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However,  the  algorithm  is  sufficiently  general  that  it  can  handle  arbitrary  nonlinear 
potential  interactions,  not  just  the  cubic  nonlinearity  of  the  Manakov  system.  We  find 
that  we  require  two  qubits  per  scalar  field  at  each  lattice  node.  The  one-particle  sec¬ 
tor  representation  of  the  Hilbert  space  allows  us  to  restrict  the  dimensionality  of  the 
wave-function  probability  amplitudes  to  just  the  dimensionality  of  the  number  of  on¬ 
site  qubits.  In  §  4  we  present  the  simulation  results  of  the  QLG  representation  for  the 
completely  integrable  Manakov  system  ( B  —  1)  as  well  as  for  the  quasi-integrable 
solitary-wave  solutions  to  the  coupled-NLS  equations  with  B  =  0.667  (this  corre¬ 
sponds  to  a  pulse  propagating  down  a  linear  birefringent  fibre).  Turbulent  soliton 
interactions  are  also  considered  using  a  potential  of  the  form  y/2/j,[\Qi\2  +  IQ2I2] 
in  the  evolution  of  the  (52-mode,  while  the  Qi-mode  satisfies  the  standard  cubic 
NLS  potential.  We  have  presented  some  preliminary  results  on  the  elastic  Manakov 
soliton  collisions  at  a  SPIE  conference  (Vahala  et  al.  2003a).  Here,  we  concentrate 
on  those  soliton  parameters  that  permit  inelastic  Manakov  soliton  collisions,  and 
consider  the  effect  of  the  cross-phase  modulation  coefficient  and  turbulence  on  the 
soliton/solitary  waves.  (We  restrict  our  use  of  the  word  ‘soliton’  to  strictly  mean 
a  form-invariant  propagating  pulse  following  nonlinear  collisions.)  Finally,  in  §  5  we 
present  our  summary  and  conclusions. 

2.  Inelastic  soliton  collisions 

(a)  Exact  2- soliton  solutions 

Consider  the  completely  integrable  Manakov  system,  equation  (1.1),  with  cross-phase 


modulation  coefficient  B  = 

1.  The  1-soliton  vector  solution  (2.4)  has  components 

Qi(x,t)  = 

f  \a  exp[-|i?  +  iTft]  sechfoR.  +  , 

(2.1) 

[i/3exp[— |.R  +  i77i]  sech[rjR  +  \R) , 

with  complex 

q  =  m.  +  iqi  =  k(x  +  i  kt). 

(2.2) 

The  parameter  R  is  real: 

„  1_[m(N2  +  I/3|2)] 

l  J’- 

(2.3) 

while  a,  (3  and  k  are  arbitrary  complex  parameters  with  £r  ^  0.  k\  determines  the 
velocity  of  the  1-soliton  components. 

Using  the  Hirota  method,  Radhakrishnan  et  al.  (1997)  have  generated  exact  2- 
soliton  vector  solutions  to  the  Manakov  system,  valid  for  all  time.  These  solutions 
are  quite  complicated.  However,  of  specific  interest  here  are  the  asymptotic  limits 
in  which  these  2-soliton  solutions  reduce  to  a  sum  of  two  widely  separated  (non¬ 
overlapping)  1-soliton  solutions 

2 

Qi(x,  t )  =  \an  exp[— |f2n  +  vq^}  sech[qnK  + 

n=l 
2 

Qz{x,  t)  =  Yl  expbi-ftn  +  i?7ni]  sech[?7nR  + 


2-^n]j 


2  Rnh 


(2.4) 
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Vn  —  VnP  4"  i-Vnl  —  kn(x  X On  4"  i&n^)j  —  hi 


KK\2  +  \Pn\2) 


with  appropriately  chosen  free  complex  parameters  an,  (5n  and  kn,  n  =  1,2.  The 
(real)  parameters  xon  predominantly  determine  the  spatial  position  of  the  peaks 
of  the  two  widely  separated  1-solitons,  kn r  the  1-soliton  amphtudes  and  kn\  the 
1-soliton  velocities. 


(6)  Asymptotic  post- collision  solitons 

The  post-collision  2-soliton  vector  solution  will  also  asymptotically  reduce  to  two 
non-overlapping  vector  1-solitons  of  the  form  of  equations  (2.4)  and  (2.5).  In  partic¬ 
ular,  for  the  choice  k\ r  >  0,  &2r  >  0,  the  post-collision  non-overlapping  1-soliton 
amplitude  ratios  are  given  by  (2.5) 

a'i  ( L  .  &i  2]  ai\(  <*iOt2  .  \  ai  2  ,  X1 

ft  U  9  Pi  J&  9Pi)\9Ptp2  {  9)  Pi  )  ’ 


where 


1  L  .  I  ^2  ,  r  * 

1~h,+  A  A+k 


g(ki,k2)  = 


-''j  (h*  —  — 
!  /  \  P*2Pl 


+  (1  +1 


h(ki,k2)  = 


fc2  4-  fcj  v  fci  +  fc|  v  ' 

Note  that  the  are  invariant  under  soliton  colhsions.  It  is  clear  that  a  certain  choice 
of  these  complex  kn  will  result  in  elastic  2-soliton  vector  collisions,  as  is  always  found 
in  scalar  soliton  collisions,  but  most  choices  of  kn  will  result  in  inelastic  2-soliton 
vector  collisions.  In  particular,  there  is  a  choice  of  parameters  for  the  kn  such  that 

a'i=0,  (2.8) 

i.e.  for  the  polarization  mode  Q\,  its  initial  non-overlapping  2-solitons  will  in  its 
post-collision  state  have  only  a  single  outgoing  soliton.  It  is  this  type  of  inelastic 
soliton  collision  that  has  encouraged  research  into  an  all-soliton  digital  information 
processor  in  a  nonlinear  optical  media  without  radiative  losses  (Jakubowski  et  al. 
1998;  Steiglitz  2001a,  b). 


3.  Quantum  lattice  gas  representation  for  coupled-NLS 

Here  we  generalize  our  QLG  algorithm  for  the  scalar  NLS  equation  (Vahala  et  al. 
2003a)  to  the  coupled  vector  NLS  system  for  arbitrary  coupling  potentials  Vj  and 
V2, 

4-  +  IQ2I]  Qi  =  o, 

4-  4-  IQii]  •  Qi  =  0. 

Our  NLS  formulation  was  based  on  the  one  particle  sector  representation  introduced 
by  Yepez  &  Boghosiah  (2002)  in  their  quantum  algorithm  for  the  (linear)  Schrodinger 
equation.  The  spatial  domain  is  discretized  into  L  spatial  nodes  and  two  qubits  per 
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scalar  wave  function  are  introduced  at  each  node  xe,  £  =  1, ...  ,L,  i.e.  qubits  x) 
for  the  scalar  field  Q i  and  qubits  |g|3)  for  Q2  at.  X(.  A  unitary  collision  operator 
that  will  recover  equation  (3.1)  is  a  factored  collision  operator  that  only  entangles 
on-site  qubits,  and  actually  only  the  on-site  qubits,  pairwise,  corresponding  to  that 
particular  scalar  field.  The  on-site  coupling  of  all  the  four  qubits  | qf),  i  =  1, . . . ,  4  is 
attained  by  introducing  a  phase  transformation  of  the  wave  functions. 

In  the  number  representation  of  the  1-particle  sector,  the  basis  set  factors  into 

lnoni  ‘ '  'n0nl)ln2n3  ' '  'n2nk)->  (3-2) 

where  only  one  of  ne0  1  is  non-zero  and  equals  unity,  and  similarly  for  n\  3.  Because 
of  the  factorization  (see  equation  (3.2)),  we  need  to  consider  the  application  of  the 
collision  and  streaming  operators  on  each  qubit  pair  separately. 

(a)  Unitary  collision  operator  for  the  qubit  pair  |go,i) 

To  evolve  the  wave  function  for  the  two  qubits  |go,i)j  a  global  quantum  unitary 
operator  C  is  constructed  from  a  tensor  product  of  quantum  gates,  each  indepen¬ 
dently  applied  to  each  spatial  site 


Ui  locally  entangles  the  two  qubits  | go,  1 )  at  each  site  and  a  local  equilibrium  can  be 
associated  with  this  on-site  unitary  collision  operator  if  \v)  is  an  eigenvector  of  Ug 
with  unit  eigenvalue:  Ui\v)  =  \v). 

The  Schrodinger  part  of  the  NLS  equation  is  recovered  by  the  VS  WAP  gate  Ue  =  U 
on  a  site-by-site  basis  (2.7) 


0  |(1  —  i)  i(l  +  i)  0 

^nls=  ?  ?  •  (3.4) 

0  i(l-fi)  |(1  -  i)  0 

\0  0  0  1/ 

with  the  Hamiltonian  representation 

Unls  =  exp[|i7r]  exp[-|i7r(c7*<7^  +  <7*0^  +  o\a2z)\  (3.5) 

written  in  terms  of  the  tensor  products  of  the  Pauli  spin  matrices 

"*=(1  0)-  <7»=(i  0).  "-“(o  -1) 

for  qubits  ‘0’  and  T’.  Note  that  br^Ls  =  the  identity  operator,  so  (7^LSIV  =  \v). 


1  0 
0  -1 


( b )  Unitary  streaming  operator  for  the  qubit  pair  |<?o,i) 

The  next  step  of  the  QLG  algorithm  is  to  stream  the  post-collision  on-site  ket 
for  qubits  |go,i)  to  nearest-neighbour  sites.  The  (unitary)  streaming  operator  Si  is 
defined  as  a  global  shift  to  the  right  of  qubit  |qo)  on  each  lattice  node,  i.e. 


Si  =  JJ  X0£-*0t+l  j 
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where  Xor-t-or+i  is  independent  of  t.  The  qubit  |<?i)  is  not  streamed  by  S\.  Thus, 
for  the  two-qubit  pair  | <70.1)1  U  entangles  the  on-site  qubit  pair,  while  x  operates  on 
the.  first,  qiibit.  of  neighbouring  sites.  Hence^  the  total  collision  matrix,  C.  does  not 
commute  with  the  full  streaming  operator  S\ . 

To  symmetrize  the  algorithm,  we  also  introduce  the  streaming  operator  S2,  which 
performs  a  global  shift  to  the  right  of  the  second  qubit  on  all  the  lattice  nodes  for 
the  qubit  pair  |<?o,i): 

L 

S 2  =  JJxit.ir+i-  (3.7) 

1= 1 

(c)  Introduction  of  the  potential  field  for  the  evolution 
of  the  wave  function  for  qubit  pair  |<3o,i) 

It  is  known  (Bialynicki-Birula  1994;  Yepez  &  Boghosian  2002)  that  the  effect  of 
an  external  potential  V\{x)  can  be  modelled  by  the  introduction  of  a  local  phase 
change  to  the  system  wave  function  for  the  qubit  pair  [go,i)-  Let  Qi  denote  the  wave 
function  for  the  qubit  pair  | Qo,  1 ) -  Thus, 

Qi(x,t)  ->  exp[iFi(a;)At]<5i(a:,t),  (3-8) 

where  At  is  the  time  advancement  after  each  step  of  the  algorithm. _ 

We  proceed  similarly  with  the  qubit  pair  (52,3),  introducing  the  \/SWAP  gate  that 
entangles  the  qubits  |g2,3)  and  the  streaming  operator  Si  that  shifts  qubit  \q2)  and 
§2  that  shifts  qubit  |<73).  One  can  now  introduce  another  external  potential  V2(x)  by 
a  local  phase  change  in  the  system  wave  function  Q2  for  the  qubit  pair  |g2,3): 

Q2(x,t)  -f  exp[iV2(x)At]Q2(x,t).  (3.9) 

(d)  Quantum  algorithm  for  coupled-NLS  equations 

The  collide-stream  sequence  of  unitary  operators  is  applied  simultaneously  to  the 
qubit  pairs  |q0,i)  and  |g2j3): 

(Qi(t  +  At)\  =  f[S%6  S2C  SJC  S2C-S?C  SiC  S??C  &6]Qi(*)\ 

\Q2(t  +  At)J  SAC  SjC  SiC-SjC  S$C  g$C  S3C}Q2(t)J  ’ 

„  •  (3-10) 

where  Sj  is  the  transpose  of  Si,  with  S?Si  =  I,  i  =  1,4  and  C  is  the  tensor  product 
of  the  unitary  collision  operator  UnlSj  equation  (3.5).  After  performing  the  collide- 
stream  sequence  (equation  (3.10)),  one  then  introduces  the  phase  change  transfor¬ 
mations  (3.8),  (3.9): 

(Qi\  ( exp[Fi(|Qi|,  |Q2|)At]  •  Qi 

Vexp[F2(|Q1|,|g2|)At]-Q2 

where  the  potential  fields  are  required  to  be  functions  of  the  wave  functions  them¬ 
selves,  Vi  =  Vi(\Qi\,  |g2|). 

The  continuum  limit  is  defined  by  scaling  the  spatial  shift  between  neigh¬ 
bouring  nodes  to  be  0{e),  the  time  advancement  At  =  0(e2)  and  the  potentials 
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Figure  1.  The  propagation  of  an  initial  (to  =  0)  exact  2-soliton  vector  solution  to  the  fully 
integrate  Manakov  system  (B  =  1)  for  the  polarization  modes  (a)  |Qi|,  ( b )  \Q2\.  The  solitons 
propagate  towards  each  other  (ti  =  25  K)  and  undergo  their  first  collision  at  t2  =  50  K. 


Figure  2.  The  post-collision  polarization  modes  for  soliton  parameters,  chosen  so  that  there  is  no 
left-propagating  soliton  for  mode  Qi  (see  equation  (2.8)).  t2+  =  51  K,  t3  =  75  K  and  U  =  100  K. 
The  intensity  in  each  mode  f  6x  \Q,(x,  t) \2  =  const,  is  preserved  to  better  than  one  part  in  10~7. 
Within  each  mode,  there  is  redistribution  of  the  intensity  in  the  outgoing  asymptotic  1-soliton 
states. 


Vi  —  s2Vi[\Qi\,  \Q2\}.  In  the  limit  e  -4  0,  it  can  be  shown  using  Mathematica  that 
the  equation  (3.10),  (3.11)  sequence  yields  the  coupled-NLS  equations 

+  +  vi[\Qi\AQ2\]  ■  Qi  =  0  +  0(e2), 

i~§tL  +  +  VzllQzl’  l^ill  •  Qi  =  0  +  0(e2), 

with  error  of  0(e2).  Equations  (3.12)  hold  for  any  choice  of  potentials  Vj  and  V2. 

It  is  also  clear  how  to  extend  this  analysis  to  a  system  of  IV-coupled-NLS  equations. 
IV-coupled  NLS  equations  arise  in  the  study  of  beam  propagation  in  a  Kerr-like 


►  as  e  ->  0,  (3.12) 

> 
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Figure  3.  The  x-t  plot  of  the  2-soliton  vector  solution  to  the  Manakov  equations  under  periodic 
boundary  conditions  and  with  the  same  parameters  as  for  the  solutions  shown  in  figures  1  and  2. 
Soliton-soliton  collisions  occur  around  t  =  50  K,  125  K  and  200  K  iterations.  There  is  a  totally 
inelastic  soliton  state  for  the  <5 1 -mode  following  the  first  collision.  Following  the  second  collision, 
the  outgoing  state  for  Qi  once  again  has  two  solitons. 

photorefractive  medium,  which  typically  exhibits  very  strong  nonlinear  effects  with 
extremely  low  optical  powers.  Partially  coherent  solitons  have  been  observed  through 
excitation  by  partially  coherent  light  (Mitchell  et  al.  1997)  as  well  as  by  an  ordinary 
incandescent  light  bulb  (Mitchell  &  Segev  1997). 

It  should  also  be  noted  that  using  the  1-particle  sector  allows  the  dimensionality 
of  the  unitary  collision  matrix  to  scale  linearly  with  the  number  of  locally  entangled 
qubits,  rather  than  scale  exponentially  with  the  total  number  of  on-site  qubits. 

4.  Quantum  lattice  gas  simulation  results 
on  inelastic  soliton  collisions 

(a)  Inelastic  Manakov  vector  soliton  collisions 

For  convenience,  we  consider  periodic  boundary  conditions  on  a  spatial  lattice  of 
6000  nodes  for  the  integrable  Manakov  equations  with  cross-phase  modulation  coef¬ 
ficient  B  =  1.  The  simulation  parameters  are  chosen  so  that  the  2-vector  solitons 
for  the  polarization  modes  Q\  and  Q 2  are  initially  non-overlapping  single  solitons, 
slightly  displaced  from  x  =  1000  and  x  —  5000,  respectively,  and  moving  towards 
each  other  for  a  collisional  interaction  (at  to  =  0  in  figure  1).  Moreover,  the  param¬ 
eters  are  chosen  so  that  the  theoretical  results  of  equations  (2.6)-(2.8)  yield  the  loss 
of  one  of  these  solitons  in  the  Qi-mode  in  the  post-collision  state. 

In  figures  1  and  2  we  show  the  results  of  the  QLG  algorithm,  and  indeed  find  the 
totally  inelastic  post-collision  state  in  which  the  second  soliton  in  mode  Q\  is  absent. 
Moreover,  in  the  fully  integrable  Manakov  system,  the  total  intensity  in  each  propa¬ 
gation  mode  Qi  is  conserved,  while  its  distribution  among  the  individual  solitons  in 
that  mode  can  be  redistributed  by  soliton  collision  (see  equations  (2.6),  (2.7)).  Our 
algorithm  conserves  this  total  individual  mode  intensity  to  better  than  one  part  in 
107,  even  in  the  single-soliton  post-collision  state.  With  periodic  boundaries,  there 
will  be  further  soliton-soliton  collisions,  as  seen  in  the  x-t  plot  in  figure  3.  Between 
the  first  and  second  soliton-soliton  collision  (52  K  <  t  <  125  K),  there  is  a  totally 
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Figure  4.  The  detailed  vector  soliton  solutions  for  modes  Qi  (dashed  curves)  and  Q?  (solid 
curves),  following  a  second  collision  for  the  parameters  given  in  figures  1-3  (£5  =  150  K, 
t6  =  175  K). 


x  x 

Figure  5.  The  propagation  of  the  same  initial  pulses  as  in  figure  1,  but  now  down  a  linearly 
birefringent  fibre  in  which  the  cross-phase  modulation  coefficient  B  =  |.  For  these  fibres,  the 
coupled-NLS  equations  are  non-integrable.  The  pulses  no  longer  propagate  as  solitons,  since 
their  amplitudes  are  no  longer  invariant  between  collisions;  rather,  they  propagate  as  soli¬ 
tary  waves  with  the  emission  of  small-scale  radiation.  Again,  the  total  intensity  in  each  mode 
Jdx  |Q<(r,t)|2  =  const,  is  preserved  to  better  than  one  part  in  10~7.  Note  that  there  is  a 
left-travelling  solitary  wave  for  the  mode  Qi  following  the  first  collision  (£4+  =  85  K),  unlike  the 
result  for  the  fully  integrable  (B  =  1)  Manakov  system. 

inelastic  collision  for  Qi.  However,  the  second  1-soliton  solution  reappears  following 
the  second  soliton-soliton  collision,  as  shown  in  detail  in  figure  4. 

( b )  Coupled  vector  NLS  of  a  linearly  birefringent  fibre 

In  a  linearly  birefringent  fibre,  the  cross-phase  modulation  coefficient  B  =  §  in 
equation  (1.1),  resulting  in  a  non-integrable  coupled-NLS  equations.  Since  there  is 
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Figure  6.  A  detailed  comparison  of  the  pulse  shapes  for  Qi  following  (a)  the  first  collision,  (6)  the 
second  collision  for  soliton  solutions  (B  —  1,  dashed  curves)  and  solitary- wave  solution  {B  =  |, 
solid  curves).  Note  the  absence  in  (a)  of  the  left-moving  soliton  for  the  fully  integrable  Manakov 
system  (B  =  1). 


k 


Figure  7.  The  intensity  spectra  for  the  Qi  mode  for  both  the  fully  integrable  Manakov  system 
(B  =  1)  at  t  =  75  K  and  t  =  150  K  and  the  mode  propagation  down  a  linearly  birefringent  fibre 
{B  —  |)  at  t  =  150  K.  The  spectrum  is  narrower  for  the  Manakov  system  ( B  =  1)  at  75  K,  since 
this  is  after  the  first  collision  and  there  is  only  one  propagating  soliton.  Following  the  second 
collision,  there  are  two  propagating  solitons  and  the  spectrum  is  similar  to  the  initial  spectrum. 
For  the  linearly  birefringent  medium  (B  =  |),  the  solitary  waves  have  a  similar  shape  to  the 
soliton  solutions  (see  figure  66)  but  the  emission  of  small-scale  radiation  gives  the  high-fc  part 
of  the  intensity  spectrum. 

only  a  slight  variation  in  B  from  B  =  1,  one  would  expect  a  quasi- integrable  solution. 
This  is  verified,  up  to  a  point,  by  our  quantum  lattice  algorithm  (figures  5  and  6). 
There  is  a  slight  decay  in  all  the  solitary- wave  amplitudes  as  they  propagate  towards 
each  other  (to  -»  ti),  with  emission  of  radiation.  In  the  post-collision  state  (t4+)  note 
that  there  is  a  quite  strong  solitary  pulse  in  the  Qi-mode  that  is  totally  absent  in 
the  fully  integrable  Manakov  system  (5  =  1).  This,  and  the  emitted  radiation  for 
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Figure  8.  A  coupled  cubic  NLS-, turbulent  NLS  system  with  <3i  integrable  (£  =  I)  except  for  its 
coupling  to  the  turbulent  <32-mode.  The  same  initial  conditions  and  parameters  as  in  figure  1 
apply.  Q i  exhibits  solitary-wave  behaviour,  including  a  substantial  pulse  at  t  =  80  K  that  would 
be  absent  in  a  Manakov  system.  The  background  radiation  is  quite  strong,  due  to  the  coupling 
to.  the  turbulent  <32-mode.  The  (^2-mode  exhibits  d-function-like  spikes  together  with  a  short 
wavelength  solitary-wave  turbulent  bath.  Collisions  between  members  of  this  bath,  as  well  as 
with  the  5-function  spikes,  are  soliton-like.  There  is  a  locking  between  the  spatial  locations  of 
the  solitary  waves  in  Qi  and  the  location  of  the  5-function  spikes  in  Qi. 


k  k 


Figure  9.  The  intensity  spectrum  for  the  Qi-  and  <52-modes  for  the  coupled  cubic  NLS-turbulent 
NLS  system.  The  solitary-wave  structure  of  \Q\  \  results  in  the  k~2  spectrum  for  medium  scales, 
while  the  background  radiation  results  in  a  k~ 5  small-scale  spectrum.  The  turbulent  |<32[-mode 
exhibits  white  noise  for  medium  scales  and  a  A:-6  small-scale  spectrum. 


the  non-integrable  case,  is  further  highlighted  in  figure  6,  where  we  compare  the 
post-collision  states  of  both  polarization  modes  for  B  =  1  and  B  =  0.667. 

Since  it  is  the  Qi  mode  that  is  the  most  affected  by  the  change  in  the  cross¬ 
phase  modulation  coefficient  B,  we  examine  how  its  intensity  spectrum  varies  in 
time  (figure  7).  First,  consider  the  fully  integrable  Manalcov  system,  B  =  1.  As 
expected,  the  spectrum  at  t  =  0  and  t  =  150  K  (which  is  after  the  second  soliton 

Phil.  Trans.  R.  Soc.  Lond.  A  (2004) 


04TA2119/12 


G.  Vahala,  L.  Vahala  and  J.  Yepez 

collision  and  the  reappearance  of  the  second  soliton  in  the  Qi-mode)  he  on  top  of 
each  other.  This  is  because  for  the  fully  integrable  Manakov  system  we  have  simply 
a  translation  in  the  position  of  the  form-invariant  solitons  which  plays  no  role  in 
the  Fourier  transform  of  \Qi\.  The  intensity  spectrum  at  t  =  75  K  (for  B  =  1)  is 
significantly  modified,  since  at  this  time  there  has  only  been  one  soliton  collision 
and  the  post-collision  state  for  |Qi  consists  of  only  one  soliton.  For  the  linearly 
birefringent  fibre  ( B  =  0.667),  we  find  solitary  pulses  and  small-scale  background 
radiation.  The  intensity  spectrum  is  statistically  independent  of  time,  and  the  effects 
of  the  radiation  emission  are  seen  in  the  raised  tail  for  wavenumbers  k  >  140.  For 
the  Q2-mode.  the  effects  of  variations  in  the  cross-phase  modulation  coefficient  B 
are  not  as  pronounced  with  a  lower  level  of  radiation'  emission.  As  a  result,  we  find 
a  raised  tail  in  the  intensity  spectrum  for  wavenumbers  k  >  275. 


(c)  Quasi-turbulent  soliton  interactions 
We  now  consider  the  following  coupled-NLS  system: 


+  M\Qi?  +  \Qz\2\  •  Qi  =  0, 

i-— ^  +  2vV\Q2\2  +  IQi I2  •  Q2  =  0. 


(4.1) 


Q 1  would  be  governed  by  the  integrable  NLS  except  for  its  coupling  to  the  Q2-mode, 
whose  evolution  is  non-integrable.  Again  we  choose  the  same  initial  profiles  (equa¬ 
tion  (2.4))  and  parameters  as  before.  For  the  mode  |Qi|,  the  initial  non-overlapping 
1-solitons  at  t  =  0  propagate  towards  each  other  like  solitary  waves  with  emission 
of  radiation.  By  t  =  80  K  (figure  8a),  the  solitary  waves  have  considerable  width 
for  this  quasi-integrable  mode  Qi-  The  evolution  of  mode  Q 2  is  considerably  dif-. 
ferent  because  of  the  underlying  turbulent-evolution  equation.  The  major  peaks  are 
extremely  narrow  (there  are  three  at  t  =  80  K;  see  figure  86),  with  the  spatial  loca¬ 
tions  of  the  two  largest  ^-function  peaks  locked  to  those  of  the  two  solitary  waves 
of  the  Qi-mode.  The  background  emission  of  Q 2  is  also  considerably  different  from 
the  Qi-mode:  it  takes  the  form  of  a  quasi-solitary  turbulent  bath  of  very  short  scales 
which  interact  with  other  structures  such  as  quasi-solitons.  This  can  be  readily  seen 
in  the  intensity  spectra  of  these  modes  (figure  9).  For  the  Qi-mode,  there  are  two  dis¬ 
tinct  spectral  regions:  a  k~2  spectrum  for  medium  range  wavenumbers,  and  a  steeper, 
k~ 5,  spectrum  for  short  spatial  scales.  The  turbulent  Q2-mode  does  not  exhibit  an 
intermediate  spectral  range,  but  it  is  fiat  (like  white  noise).  For  very  short  scales  we 
see  an  intensity  spectral  decay  of  k~6. 


5.  Conclusion 

A  lattice-based  quantum  representation  is  developed  for  the  solution  of  coupled- 
NLS  equations:  equations  that  can  model  the  propagation  of  pulses  down  a  birefrin¬ 
gent  fibre.  The  algorithm  requires  two  qubits  per  spatial  node  for  each  polarization 
mode  and  introduces  quantum  entanglement  on  the  on-site  qubits.  This  entangle¬ 
ment  becomes  strongly  non-local  due  to  the  application  of  the  streaming  operator.  A 
symmetrized  form  of  streaming  on  each  qubits  results  in  a  more  accurate  numerical 
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scheme  (equation  (3.12)),  and  in  finite-difference  form  leads  to  an  unconditionally  sta¬ 
ble  algorithm.  The  nonlinear  interaction  between  the  polarization  modes  is  achieved 
by  appropriate  phase  transformations  on  all  four  on-site  qubits,  and  we  have  tested 
our  algorithm  on  a  particularly  exacting  problem  of  totally  inelastic  vector  soliton 
collisions  of  the  completely  integrable  Manakov  system. 

In  this  vector  NLS  algorithm,  quantum  phase  coherence  needs  to  be  enforced 
for  only  16  quantum-gate  operations,  after  which  measurements  are  required  so 
as  to  determine  the  amplitudes  of  the  wave  functions.  These  amplitudes  are  then 
used  to  construct  the  nonlinear  interaction  potential.  It  is  clear  that  our  quantum 
algorithm  for  vector  solitons  would  be  difficult  to  implement  experimentally.  Even 
though  it  is  a  one-dimensional  model,  it  is  complex  because  amplitudes  must  be 
extracted  from  the  quantum  measurement  process  and  not  simply  from  probabili¬ 
ties  which  are  more  naturally  extracted  from  the  projective  von  Neumann  quantum 
measurement  process  over  a  large  ensemble  of  molecules,  as  occurs  in  nuclear  mag¬ 
netic  resonance  spectroscopy  (Pravia  et  al.  2002,  2003).  Nevertheless,  it  is  hoped 
that  this  one-dimensional  quantum  algorithm  will  serve  as  a  stepping  stone  towards 
more  complicated  quantum  algorithms  with  sufficiently  many  qubits  per  site  that 
the  experimental  trade-off  between  the  difficulty  of  achieving  quantum  entanglement 
and  the  limitation  of  measurable  classical  information  that  can  be  extracted  from 
the  quantum  computer  begins  to  work  in  our  favour.  Finally,  it  is  interesting  to  note 
that  the  QLG  representation  is  very  efficiently  parallelized  on  a  classical  computer. 

This  work  was  supported  by  the  Directorate  of  Computational  Mathematics,  US  Air  Force  Office 
of  Scientific  Research. 
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